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Abstract 

We propose a two-dimensional model consisting of photons and photon pairs. 
In the model, the mixed gas of photons and photon pairs is formally equiva- 
lent to a two-dimensional system of massive bosons with non- vanishing chem- 
ical potential, which implies the existence of two possible condensate phases. 
Using the variational method, we discuss the quantum phase transition of 
the mixed gas and obtain the critical coupling line analytically. Moreover, 
we also find that the phase transition of the photon gas can be interpreted 
as second harmonic generation. We then discuss the entanglement between 
photons and photon pairs. Additionally, we also illustrate how the entangle- 
ment between photons and photon pairs can be associated with the phase 
transition of the system. 

Keywords: Bose-Einstein condensate, Quantum phase transition, 
Entanglement 

1. Introduction 

Bose-Einstein condensate (BEC) is the remarkable state of matter that 
spontaneously emerges when a system of bosons becomes cold enough that 
a significant fraction of them condenses into a single quantum state to min- 
imize the system's free energy. Particles in that state then act collectively 
as a coherent wave. The phase transition for an atomic gas was first pre- 
dicted by Einstein in 1924 and experimentally confirmed with the discovery 
of superfluid helium-4 in 1938. Obviously, atoms aren't the only option for a 
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BEC. In recent years, with the development of techniques, the phenomenon 
of BEC was observed in several physical system [1-9], including exciton po- 
laritions, solid-state quasiparticles and so on. We know that photons are the 
simplest of bosons, so that it would seem that they could in principle un- 
dergo this kind of condensation. The difficulty is that in the usual blackbody 
configuration, which consists of an empty three-dimensional (3D) cavity, the 
photon is massless and its chemical potential is zero, so that the BEC of 
photons under these circumstances would seem to be impossible. However, 
very recently, J. Klaers, etc. have overcome both obstacles using a simple 
approach [10,11]: By confining laser light within a two-dimensional (2D) cav- 
ity bounded by two concave mirrors, they create the conditions required for 
light to thermally equilibrate as a gas of conserved particles rather than as 
ordinary blackbody radiation. 

What is more, it is well known that there are many fascinated optical 
effects in the nonlinear medium, for instance, reduced fluctuation in one 
quadrature (squeezing) [12], sub-Poissonian statistics of the radiation field 
[13], or the collapse-revivals phenomenon [14]. Especially, in the nonlinear 
medium, a photon from the laser beam can couple with other photons to 
form a photon-pair (PP) [15-18]. The essence of PP has been investigated 
by many authors [19-21]. However, inspired by the experimental discovery 
of BEC of photons, in this letter we construct another interesting 2D model 
consisting of photons and PPs. In this model, the mixed system of photons 
and PPs is formally equivalent to a 2D gas of massive bosons with non- 
vanishing chemical potential, which implies the existence of two possible 
condensate phases, the mixed photon-PP condensate phase and the pure 
PP condensate phase. By means of a variational method we investigate the 
quantum phase transition of the mixed photon gas. Especially, we find that 
the quantum phase transition of the photon gas can be interpreted as second 
harmonic generation. We then discuss the entanglement between photons 
and PPs. By investigating the entanglement in the ground state and the 
dynamics of entanglement, we also illustrate how the entanglement between 
photons and PPs can be associated with the phase transition of the system. 
The investigation of these questions is important both for its connection with 
quantum optics and for its practical applications to harmonic generation and 
quantum information. 

The remainder of this paper is organized as follows: In Sec. II, we theoret- 
ically investigate the phenomenon of BEC of photons and PPs in a 2D optical 
microcavity. The entanglement between photons and PPs is investigated in 
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Sec. III. Finally, we make a simple conclusion. 

2. Bose-Einstein condensation of photons and photon pairs 

2.1. Description of the photon pair 

We start with the description of the PP. History speaking, the essence 
of the PP is presently still under discussion [19-21], and there exist many 
different ways to obtain it. However, here we will use the standard procedure 
[22] in the construction of harmonic generation to derive the PP. We know 
that the presence of an electromagnetic field in the nonlinear material causes 
a polarization of the medium and the polarization can be expanded in powers 
of the instantaneous electric field: 

P(r, t) = x (1) E(r, t) + X (2) E 2 (r, t) + .... (1) 

Here, the first term defines the usual linear susceptibility, and the second 
term defines the lowest order nonlinear susceptibility. Ignoring the high order 
parts (i.e. only expend the polarization to second order in electric field E), 
we find that the Hamiltonian describing the interaction of the radiation field 
with the dielectric medium is decomposed into two terms: 

H m t Hn ne -\- H non n ne 
H lme = -jx^E 2 (r,t)dr 

Hnonline = - j X ^E 3 (r,t)dr 

where Hi ine represents the energy of the linear interaction and H non [ ine the 
nonlinear interaction. 

It is well known that the electric field operator in a microcavity can be 
expanded in terms of normal modes [23] as 

/ ft, > \ 1/2 

E(r, = (^j (a k e-^ k ' r + a+e^^) , (3) 

where and are the annihilation and creation operators of photons with 
frequency Ufa, and they all obey the usual boson commutation rules. V is the 
normalization volume, e is the dielectric constant of the medium and e k is 
the unit polarization vector with the usual polarization indices omitted for 
simplicity. Substituting (3) into (2), for the linear interaction part, we find 
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(2) 



that it consists of two processes, dissipation and two-photon absorption(or 
emission). Here, dissipation is essentially also a two-photon process, in which 
one photon is absorbed by the medium, meanwhile another one is emitted. 
The linear interaction can be ignored, if the incident photon field frequency 
uo is well below the electronic transition frequencies of the medium. In that 
case, we need only consider the nonlinear interaction, which has the simple 
form 

H nonline = —j= Xk,k' ( &k+k'^kOk' + H.C.\ (4) 

under the requirements of phase matching. Above, the operator a represents 
the normal photons, b represents the coupling PP, and where Xk,k' is the 
coupling matrix element. The interaction energy in (4) consists of two terms. 
The first term &£ +k /akOk' describes the process in which two normal photon 
with wave- vector k and k' couple into a PP with wave- vector K = k + k', 
and the second term describe the opposite process. The energy is conserved 
in both the processes. 

2.2. Free-photon dispersion relation inside the optical microcavity 

In this letter, we restrict out investigation inside a 2D optical microcavity. 
The microcavity, as shown in Fig. 1, consists of two curved dielectric mirrors 
with high reflectivity (about 99.9), which ensure prefect reflection of the 
longitudinal component of the electromagnetic field within the cavity. In 
addition, the transverse size of the cavity is much larger than its longitudinal 
one. 

We know that for a free photon, its frequency as a function of transversal 
(k r ) and longitudinal (k z ) wave number is co = c [k% + ^] 1//2 . However, in the 
case of photons confined inside the microcavity, the vanishing of the electric 
field at the reflecting surfaces of the curved-mirrors imposes a quantization 
condition on the longitudinal mode number k z , k z = nir/D(r), where n is 
an integer and where D(r) = D$ — 2(R — \JR 2 — r 2 ) is the separation of 
two curved-mirrors at distance r from the optical axis, with Dq the mirror 
separation at distance r = and R the radius of curvature. 

In the present work, we consider to fix the longitudinal mode number of 
photons by inserting a circular filter into the cavity. The filter is filled with 
a dye solution, in which photons are repeatedly absorbed and re-emitted 
by the dye molecules. Thus, it also plays the role of photon reservoir. We 
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know that the longitudinal size of the cavity (i.e., the distance between the 
mirrors) is very small. The small distance D(r) between the mirrors causes 
a large frequency spacing between adjacent longitudinal modes, comparable 
with the spectral width of the dye. Modify spontaneous emission such that 
the emission of photons with a given longitudinal mode number, n = q in our 
case, dominates over other emission processes. In this way, the longitudinal 
mode number is frozen out. For fixed longitudinal mode number q and in 
paraxial approximation (r <C R, k r <C k z ), we also find that the dispersion 
relation of photons approximatively becomes u ~ qiic/D + ck 2 D /2qii. The 
above frequency-wavevector relation, upon multiplication by h, becomes the 
energy-momentum relation for the photon 



E » m ph c 2 + ^L, (5) 
2m ph 

where m p h = hq7r/D c = fvuj e ff/c 2 is the effective mass of the confined pho- 
tons. At low temperatures, it is convenient to redefine the zero of energy, so 
that only the effective kinetic energy, 

E«%£, (6) 
2m ph 

remains. The above analysis shows that for the photon confined inside the 2D 
microcavity, it is formally equivalent to a general boson having an effective 
mass m p h = hu e ff/c 2 , that is moving in the transverse resonator plane. 

Furthermore, we here consider the case that the microcavity (except the 
filter part) is filled with a Kerr nonlinear medium exhibiting significant third- 
order optical nonlinearity. Due to the nonlinear effect, photons can couple 
into PPs. If we connect the non-vanishing effective photon mass to the 
previous analysis of the PPs, in this case we then can rewrite the nonlinear 
interaction H nonHne as 



nonline 



-7= X^Vr ( 6 k r +kv a k r akv + H.c) (7) 

where S is the surface area of the 2D cavity, and where a-^ r and a k / r are 
the annihilation operators of massive photons with transverse wavevectors 
k r and k r /, respectively, and &k r +k' r are ^ ne creation operator of the massive 
PPs with transverse wave- vector K r = k r + k' r . Here, it should be remarked 
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that the existence of effective photon mass makes the thermodynamics of 
this 2D mixed gas of photons and PPs different from the usual 3D photon 
gas. For the 2D system, thermalization is achieved in a photon-number- 
conserving way (N = N a + 2N b ) with nonvanishing chemical potential fi, 
by multiple scattering with the dye molecules, which acts as heat bath and 
equilibrates the transverse modal degrees of freedom of the photon gas to the 
temperature of dye molecules. 

2.3. BEC of photons and photon pairs 

In virtue of the above analysis, we consider the following basic Hamilto- 
nian, to give a simple model of PP formation (with h — 1 throughout this 
letter) 

H fl = H a + H b + H ab , (8) 

with 

H a = E £: a k a k + Y E 4+k'-k» a k» a k'ak 

k p k,k',k" 

H b = E(^r h -^)KK + ] f E K+u-uft'bvbK (9) 

k v 7 k,k',k" 

H a b = ^Yl a£K> b wa k - -j- E Xk,k'( & k+k' a k'«k + H.c.) 
k,k' v k,k' 

Above, H a and H b denote the pure photon and PP contributions, and H ab 
refers to the interaction between them. In the dilute gas limit u a>ai u bjb , u ajb 
are proportional to the two-body s-wave photon-photon, photon-PP, and PP- 
PP scattering lengths [23], respectively, and Xk,k' characterizes the coupling 
strength, encoding that PPs are composed of two massive photons. Note that 
in (9) we ignore the chemical potential term fiN a from the Hamiltonian [24] . 
Additionally, in (9), we also drop the subscript of the transverse wave- vectors 
of photons and PPs, i.e. we take k r = k. 

It is well known that for a general massive Boson-Boson pairs gas, at 
absolute zero temperature, there exists a BEC consisting of two possible 
condensate phases [25-27]: (i) Both the single boson and the pair of bosons 
are condensed, (ii) The pair of bosons are condensed but the single boson is 
not. Now we know that in the case of photons (or PPs) confined inside the 
microcavity, the subsystem of photons (or PPs) is formally equivalent to a 
2D gas of massive bosons with non-vanishing chemical potential. Thus, the 
feature should also survive for the 2D mixed gas of photons and PPs. In 
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the condensate phase, a macroscopic number of particles occupy the zero- 
momentum state, and it is useful to separate out the condensate modes from 
the Hamiltonian. Follow the process, we find that at the BEC state, the 
grand canonical Hamiltonian of the mixed system has the form 

H^ = H + 5H, (10) 

where 

H = -/i 6 6+6 + Y (« + ^) 2 + {b+bf 
+ ^ a + a b+b- ^(a+a+b + b+aa) 

is the condensate part of the Hamiltonian with /i b = 2fi + u bb /S the modified 
chemical potential of PPs, and where 5H = H (ak, 6k) is the perturbation 
part and it is a complex function of the non-condensate modes. Here, we 
mention that at the condensate state we only need consider the case of sin- 
gle mode coupling, thus we can treat the coupling matrix element x as an 
adjustable constant. Note that in (11) we also drop the zero- momentum 
subscript of the creation and annihilation operators of photons and PPs. In 
the present work, we mainly aim to investigate the phenomenon of BEC of 
photons and PPs, thus, hereafter we will ignore the perturbation part and 
approximately write the Hamiltonian as the form Hq. The Hamilto- 

nian commutes with the total photon number N = a + a + 26 + 6 and n = N/S 
is the particle density. 

Up to now we have not made a careful distinction between the two pos- 
sible condensate phases. However, for the mixed system, working out the 
ground-state phase diagram is very important. Here, we intend to employ 
the variational principle method for finding the ground-state configurations 
of the present system and examining their dependence from the microscopic 
parameters. In other word, we aim to work out the ground-state phase dia- 
gram of the mixed system starting from the study the semiclassical equation. 
Clearly, we know that the chemical potential \i of the mixed system allows the 
total photon number (whether free or bound into PPs) to fluctuate around 
some constant average value N, then the total number of photons need only 
be conserved on the average value. For convenience, hereafter we assume 
that N is an even number and thus M = N/2 denotes the maximum number 
of the PPs. Furthermore, in this case we also introduce a new operation, 
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namely the double photon creation operation with the relation 



(^ + ri0)-^(a + ) 2m |0) (12) 
c+c |0) = 2a+a\0) 

where |0) is the vacuum state. Using the new operator, we construct the 
Gross-Pitaevskii (GP) states [28] 

|^ GP > = - 7 L=[ac + + /% + ] M |0> (13) 

as the trial macroscopic state. Here, a = \a\e l0a and (3 = \(3\e t9b are com- 
plex amplitudes with \a\ 2 = N a /N and \f3\ 2 = 2N b /N the photon and PP 
densities, respectively. 9 a and 9 b (real valued) denoted the phases of each 
species. Obviously, the parameters a and (3 satisfy the normalized condition 
M 2 + \P\ 2 — 1- With the help of the GP state, the semiclassical model 
Hamiltonian H(a, (3) is given by 

H = lim ^o P \H,\p P ) 

= [-ti b \(3\ 2 + 2u aa n |a| 4 + u bh n |/3| 4 /2 , (14) 

2 i n\ 21 / fn n I 1 2 l\ n\2 



+u a bn\a\ \(3\ ] /xv / 2n — 2 |a| y\/3\ cos9 

where 9 = 9 b — 29 a is the phase difference. Considering the conserved condi- 
tion |a| 2 + |/3| 2 = 1, we next introduce a new variables s = \a\ 2 — \f3\ 2 . Using 
the new notation, we rewrite the model Hamiltonian as 

H = -As 2 - 2 7 s + f - y/2 (1 - s)(l + s) cos 9, (15) 

with 

\ _ V2n ( U ab Uaa U bb \ 

_ V2n I Ubb _ Uaa _ t±b\ 

I x v 8 2 An) 

t _ / Uaa I «o6 I "66 _ Mb 

x^2~ r 4~ r 8 n / 

According to the variational principle, we minimize the energy H(a, f3) with s 
and 9 as variational parameters. We then obtain the optimum values [i.e.(s, 
9)] of parameters for the ground state as follows: 

^ \ (s, Oom), 7 -A + l>0 ' 1 ° j 
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where — 1 < s < 1 is the solution of the equation As+7 = (3s — 1) /2a/2(1 — s) 
(the explicit value can be obtain by graphical solution method, it is gener- 
ally too messy to be shown here). The result, together with the fact that the 
parameters \a\ 2 and |/3| 2 denote the photon and PP densities, indicates that 
when 7 — A + 1 < the system converts from the mixed photon-PP phase 
to the pure PP phase. We therefore can interpret this line 7 — A + l = 0as 
the threshold coupling for the formation of a predominantly PP state. Here, 
it is to be mentioned that for the pure PP phase, s — —1 thus the relative 
phase 9 cannot be defined. 



3. Entanglement between photons and photon pairs 

3.1. Entanglement of the ground state 

At the BEC state, one may consider the mixed gas of photons and PPs as 
a bipartite system of two modes. For the present system, the entanglement 
of two modes is always closely associated with the phase transition of the 
system. Moreover, the two modes, be they spatially separated, and differing 
in some internal quantum number, are clearly distinguishable subsystems. 
Thus, the state of each mode can be characterized by its occupation number. 
By using the fact that the total number of photons N is constant, a general 
state of the system (in the Heisenberg picture) can be written for even N in 
terms of the Fock states by 

M 

\i)) = ^c m \2m,M-m), (17) 

m=0 

where m is the half population of particles in photon mode a, and c m is 
the coefficients of the state. In the Fock representation, the GP state also 

M 

can be reexpressed as \ipcp) = Yl 9m |2m, M — m), with coefficients g m = 

m=0 

a m /3 M ~ m . The standard measure of entanglement of the bipartite 



m\(M-rn)\ 

system is the entropy of entanglement S(p) 



M 



S(p) = -$> m | 2 log 2 (|c m | 2 ), (18) 



m=0 



which is the von Neumann entropy of the reduced density operator of either 
of the subsystems [29]. In the present system, the maximal entanglement also 
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can be obtain by optimizing the expression (18) with respect to \c m \ . By 

M 

imposing the normalization condition Yl \ c m\ = 1, we finally get S^ax = 

m=0 

log 2 (M + 1), which is related to the dimension M + 1 of the Hilbert space 
of the individual modes. 

Using expression (18) and the coefficients c m obtained through exact di- 
agonalization of the Hamiltonian (11) as done in the atom- molecule model 
[30], we plot in Fig. 2 the entropy of entanglement of the ground state as a 
function of the parameters A and 7. We note that in this letter we restrict 
our attentions to the repulsive case, i.e., we restrict A < throughout the 
letter. From Fig. 2, we observe that the entanglement entropy exhibits a 
sudden decrease close 7 — A + 1 = 0. This is indicative of the fact that across 
the line 7 — A + l = 0a quantum phase transition occurs. 

To gain more information associated with the quantum phase transition 
of the present system, we also depict in Fig. 3 the entropy of entanglement 
(solid line) and the expectation value (dashed line) of the scaled PP number 
operator of the ground state as a function 7 for fixed parameter value A. 
From Fig. 3, we see that the average value of the number of PPs increases 
as 7 increases. Especially, when 7 — A + 1 < 0, the average number of PPs is 
maximal. The result confirms that there indeed exists a phase transition for 
the present system in the ground state. Furthermore, we also find that the 
ground-state entanglement entropy is not maximal at the critical line, i.e. in 
the region 7 — A + 1 > 0, the system is always strongly entangled. Ref. (28) 
gives the property responsible for the long-range correlation. Additionally, 
we also consider that the trait is associated with the symmetry-broking of 
the coupling term of the system. Due to the asymmetric form of the cou- 
pling term, the pure photon condensation will be forbidden. As a result, in 
the mixed condensate phase the imbalance (2Ni — N a ) /N between the two 
modes is always very small, which is responsible for the strongly entangle- 
ment. 

In addition, if we connect the non-vanishing photon mass m p h to the 
longitudinal wave number k z by the relation m p h = hk z /c with k z = qn/DQ, 
then we find that the quantum phase transition of the photon system can be 
interpreted as second harmonic generation. When 7 — A + 1 < 0, almost all 
photons with frequency u = ck z couple into PPs with frequency uj = 2ck z . 
In this case, the entanglement between the photons and PPs is very small, 
and the entropy of entanglement is close to zero. 
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3.2. Dynamics of entanglement 

In the above analysis, we have investigated the entanglement of the 
ground state. We found that in the ground state, across the phase tran- 
sition line the entanglement entropy exhibits a sudden change. To gain a 
better understanding of the influence of ground-state phase transition to en- 
tanglement, in this subsection we investigate the dynamics of entanglement. 
In studying the dynamics of the system, we first need express a general state 
in the form of temporal evolution (i.e., need change the expression of a gen- 
eral state from Heisenberg picture to Schrodinger picture). Following the 
standard procedure, we can obtain 

|^(f)> = tf(f)hMo)> 

M (19) 
= £ c m (t)\2m,M-m) ' 1 ; 

m=0 

M 

where, U (t) — Yl IVO (V'nl ex P ( — iE n t) is the temporal operator with \ip n ) 

n=0 

the eigenstates of the system having energy E n , and |-0 (0)) is the initial 
state. Here, the time dependence of coefficients c m (t) are given by c m (t) = 
(2m, M — to | U (t) (0)). Subsequently, the entanglement entropy given in 

M 

(18) can be rewritten as S(p) = — £ l c m(0l log 2 (|c m (t)| ). In this case, 

m=0 

the entanglement entropy depends on both the choice of initial states and 
the value of microscopic parameters {A, 7}. At the present work, we consider 
that the mixed system is in the BEC state, thus here choosing GP state as 
the initial state is suitable. By adjusting the GP coefficients {\a\ 2 , \(3\ 2 } and 
the microscopic parameters {A, 7}, in this subsection we also want to know 
if the ground-state phase transition also characterizes different dynamics. 

The time evolution of the entanglement entropy for different initial state 
and interaction parameters is shown in Fig. 4. From Fig. 4 we observe the 
features of quantum dynamics, such as the collapse and revival of oscillations 
and non-periodic oscillations. Additionally, we also find that the amplitude 
of the entanglement entropy is smaller in the region 7 — A + 1 < con- 
trast with in the region 7 — A + 1 > 0. Especially, we note that the greater 
the imbalance \(3\ 2 — \a\ 2 between the two modes in the initial state, the 
clearer the difference can be observed. To understand the physical reason 
for the above phenomenon, we also need rewrite the general state given in 

M 

(19) in terms of the eigenstates of the system \ip (t)) = £ c(n,t) \ip n ), where 

n=0 
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\c(n,t)\ 2 = |c(n,0)| 2 = \(ip n \ exp(— iE n t) (0))| 2 can be explained as the 
transition probability of the system from the initial state \ip (0)) to the cor- 
responding energy eigenstates \ip n ) at any time t. We have already known 
that for the ground state across the phase transition line 7 — A + 1 = the 
entanglement entropy exhibits a sudden decrease. This is why in the region 
7 — A + 1 < the amplitude of the entanglement entropy becomes smaller. 
In addition, in this phase transition region we also investigative the depen- 
dence relation between the ground-state transition probability |c(0,t)| 2 and 
the imbalance |/3| 2 — \a\ 2 of the initial state. The result is shown in Fig. 5. 
From Fig. 5, it is obvious that with the increasing of the initial-state im- 
balance the ground-state transition probability becomes greater. Thus the 
greater imbalance between the two modes in the initial state can lead the 
clearer difference of the amplitude for different region. 

4. Conclusion 

In this work, we have proposed a 2D model consisting of photons and PPs. 
In the model, the mixed gas of photons and PPs is formally equivalent to a 
2D system of massive bosons with non-vanishing chemical potential, which 
implies the existence of two possible condensate phase. Based on the GP 
state and using the variational method, we have also discussed the quantum 
phase transition of the mixed gas and have obtained the critical coupling 
line analytically. Especially, we have found that the phase transition of the 
photon gas can be interpreted as second harmonic generation. Moreover, by 
investigating the entanglement entropy in the ground state and general state, 
we have illustrated how the entanglement between photons and PPs can be 
associated with the phase transition of the system. 
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Figure 1: Scheme of the optical microcavity: The microcavity consists of two curved- 
mirrors with high reflectivity. A filter filled with a dye solution is inserted into it, in which 
photons are repeatedly absorbed and re-emitted by the dye molecules. The other part of 
the cavity is filled with a Kerr-like nonlinear medium. 




Figure 2: Variation in the entropy of entanglement of the ground state with respect to the 
parameters A and 7. Here, we set N = 200 and u aa — u^/A — 0.25. 
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Figure 3: The average photon pair occupation number and the entanglement entropy for 
the ground state as a function of the 7 for fixed parameter value A = 0. Here, we set 
N = 200 and u aa = u bb /A = 0.25. 
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Figure 4: Time evolution of the entanglement entropy for different initial states 
tp (\a\ 2 , |/3| 2 ^ and microscopic parameters {A, 7}. Form top to bottom the GP coef- 
ficients used arc j|a| 2 ,|/3| 2 } = {0.5,0.5}, {|a| 2 ,|/3| 2 } = {0.25,0.75} and {|a| 2 ,|/3| 2 } = 

{0,1}. From left to right the microscopic parameters used are {7, A} = {—2,0} and 
{7, A} = {—0.5, 0}, respectively. Here, we set N = 20 and u aa = u^/A = 0.25. 
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Figure 5: Variation of the ground-state transition probability |c (0, t)\ 2 with respect to the 
initial-state imbalance \(3\ 2 — \a\ 2 . Here, we set N = 20 and u aa = Ubb/4 = 0.25. 
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